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SECTION  1 


INTRODUCTION 

Naturally  occurring  sources  of  inffasound  include  (but  are  not  limited  to)  severe  weather, 
volcanoes,  bolides,  earthquakes,  surf,  mountain  waves,  and,  the  focus  of  this  research,  nonlinear 
ocean  wave  interactions.  Man-made  sources  of  inffasound  also  exist,  such  as  airplane  activity, 
military  testing,  rocket  launches,  and  nuclear  explosions  (Bedard  and  Georges,  2000).  Due  to  its 
low  frequency,  inffasound  can  travel  global  distances  with  relatively  low  attenuation  while 
higher  ffequency  audible  sound  is  usually  dissipated  at  shorter  ranges  (Drob  et  al.,  2003). 
Infrasonic  waves  may  travel  through  the  atmosphere  between  the  Earth’s  surface  and  the 
thermosphere.  The  variability  of  wind  and  temperature  with  height  determines  the  advection  and 
refraction  of  inffasonic  waves  in  the  atmosphere  (Gossard  and  Hooke,  1975). 

Infrasound  station  I59US,  located  at  latitude  19.5915  N,  longitude  155.8936  W  (near  Kailua- 
Kona,  Hawaii)  and  operated  by  the  University  of  Hawaii,  is  part  of  the  global  infrasound 
network  of  the  IMS  (Fig.  1,  Vivas-Veloso  et  al.,  2002).  Inffasound  arrivals  at  I59US  (Fig.  2)  have 
been  studied  in  detail  to  characterize  the  coherent  background  sound  field  (Garces  and  Hetzer, 
2001,  2002,  2003).  A  large  portion  of  this  ambient  field  is  related  to  a  pervasive  signal  known  as 
microbaroms.  Microbaroms  are  observed  as  a  continuous  atmospheric  pressure  oscillation  with 
most  of  its  energy  between  0. 1  and  0.5  Hz  (Fig.  3).  Microbaroms  are  detected  at  the  I59US  array 
as  coherent  bursts  with  durations  of  minutes,  an  RMS  amplitude  varying  between  ~10  mPa  and 
-100  mPa,  and  apparent  phase  speeds  with  roughly  the  speed  of  sound  (-340  m/s)  (Figs.  4,5). 
For  inffasonic  stations  near  the  ocean,  microbaroms  comprise  the  low-wind  noise  floor.  The 
microbarom  peak  is  in  the  midst  of  the  detection  region  for  1-kiloton  nuclear  explosion  tests. 
Therefore,  microbaroms  can  obscure  a  signal  of  interest. 

The  objectives  of  this  project  were  to  characterize  microbarom  signals  observed  in  the  Pacific 
and  model  the  source  processes  that  generate  these  signals,  with  the  aim  of  determining 
infrasonic  detection  thresholds  in  the  microbarom  ffequency  range.  We  addressed  these 
objectives  by  studying  microbarom  signals  from  selected  case  studies  and  associating  the 
observed  infrasonic  amplitudes,  signal  arrival  direction,  and  microbarom  frequencies  with  the 
theoretical  infrasonic  source  field  generated  by  a  prescribed  ocean  surface  wave  spectrum.  In  this 
report  we  describe  how  the  amount  of  coherent  and  incoherent  infrasonic  energy  in  the 
microbarom  frequency  band  is  determined  by  a  complex  combination  of  lower  and  upper 
atmosphere  winds  and  the  proximity,  dynamics,  and  intensity  of  open-ocean  swells.  This  report 
summarizes  the  results  of  various  meeting  presentations,  a  Masters  thesis  (Willis,  2004),  and 
three  peer-reviewed  papers  (Willis  et  al.,  2004,  Garces  et  al.  2004,  Willis  et  al.,  submitted). 

Microbaroms  were  first  reported  by  Benioff  and  Gutenberg  (1939)  on  an  electromagnetic 
microbarograph  at  the  Seismological  Laboratory  of  the  California  Institute  of  Technology.  At  the 
time  of  their  studies,  there  was  no  accepted  hypothesis  for  microbaroms  or  their  better-known 
seismic  counterparts,  microseisms.  Longuet  Higgins  (1950)  was  the  first  to  recognize  and 
develop  a  theory  for  the  excitation  of  microseisms  through  the  interaction  of  standing  ocean 
waves  with  the  sea  floor.  Posmentier  (1967)  presented  a  theory  similar  to  the  approach  of 
Longuet-Higgins  that  explained  the  source  generation  of  microbaroms.  Posmentier’ s  theory 
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described  a  nonlinear  pressure  perturbation  that  arises  at  the  air-sea  interface  when  two  ocean 
waves  of  opposite  direction  and  similar  frequencies  meet.  The  corresponding  acoustic  wave  was 
shown  to  gain  properties  of  the  interfering  ocean  wave  train  where  acoustic  amplitude  is 
proportional  to  the  product  of  the  opposing  ocean  waves  and  frequency  is  twice  that  of  the 
individual  ocean  waves.  Arendt  and  Fritts  (2000)  extended  this  theory  to  an  arbitrary  spectrum  of 
ocean  surface  waves.  They  found  that  the  frequency-doubling  nonlinear  interactions  of  pairs  of 
ocean  waves  traveling  in  nearly  opposite  directions  produces  propagating  acoustic  waves  of  an 
isotropic  nature.  All  terrestrial  ocean  surface  waves  contain  phase  speeds  and  wavelengths  much 
smaller  than  acoustic  phase  speeds  and  wavelengths.  Thus,  single  ocean  waves  cannot  couple 
sound  into  the  atmosphere.  However,  sound  can  be  coupled  into  the  atmosphere  when  ocean 
waves  interact  nonlinearly  and  the  sum  of  the  horizontal  wavenumbers  of  the  ocean  surface 
waves  is  nearly  zero.  This  occurs  only  when  an  ocean  wave  spectrum  contains  components  of 
nearly  identical  frequencies  traveling  in  nearly  opposite  directions.  Ocean  waves  propagating 
with  identical  frequencies  in  opposite  directions  are  predicted  to  radiate  sound  vertically,  leading 
to  efficient  ensonification  into  the  ocean  and  seabed  but  ensuring  that  atmospheric  sound  never 
reaches  the  ground  again.  Thus,  microseisms  propagate  through  the  ground  as  a  result  of  vertical 
excitation  through  the  ocean;  microbaroms  propagate  to  infrasonic  stations  after  near  horizontal 
propagation  through  the  atmosphere  (Tebulevich,  1995). 

Studies  by  Saxer  (1945,  1954),  Daniels  (1953,  1962),  Donn  and  Posmentier  (1967),  Donn  and 
Naini  (1973)  and  Rind  (1980)  confirmed  that  the  microbarom  and  microseism  source  is  related  to 
strong  storms  over  the  ocean  and  the  resulting  high  seas.  In  addition  to  major  weather  systems 
(cold  fronts,  high  and  low  pressure  areas)  and  significant  wave  heights,  Rind  (1980)  compared 
expected  source  locations  of  microbaroms  and  microseisms  recorded  at  Palisades,  New  York 
with  dominant  wave  period  and  mean  propagation  charts  on  a  5°  grid  provided  by  the  Navy  Fleet 
Numerical  Weather  Center.  Since  microbaroms  theoretically  contain  frequencies  twice  those  of 
the  producing  ocean  waves.  Rind  attempted  to  correlate  the  period  of  the  observed  microbaroms 
with  ocean  wave  half  periods  in  the  expected  source  regions.  The  mean  propagation  charts  were 
used  to  correlate  the  expected  source  regions  with  areas  that  contained  opposing  wave  trains. 
The  methods  of  Rind  (1980)  were  not  effective  in  the  case  studies  presented  in  this  study.  It 
should  be  mentioned  that  Rind  discussed  the  error  potential  in  correlating  microbaroms  with 
mean  wave  parameters  on  such  a  coarse  grid  instead  of  using  any  entire  spectrum  of  waves. 
Advancements  have  been  made  in  ocean  wave  modeling  since  1980,  which  has  made  this  present 
research  possible. 

Strong  microbarom  signals  tend  to  arrive  from  regions  of  marine  storminess  because  marine 
storms  can  produce  regions  of  high  amplitude  ocean  waves  converging  from  opposite  directions. 
Consider  a  symmetrical,  westward  moving  tropical  cyclone  of  the  northern  hemisphere  (Fig.  6a). 
The  surface  low  will  contain  similar  fetch  characteristics  on  all  sides  of  the  center,  with  the  north 
side  of  a  westward  moving  storm  in  the  northern  hemisphere  generally  containing  the  highest 
winds  and  seas.  After  the  westward  moving  storm  passes  a  given  point,  the  swell  spectrum  at 
that  point  will  eventually  become  confused  with  wave  trains  traveling  in  opposite  directions  that 
were  created  on  both  the  northeast  and  southwest  sides  of  the  low  center  (Fig.  6b).  Thus,  the 
theoretical  guidelines  for  microbarom  productions  will  be  met  in  the  storm  wake  region,  which  is 
east  of  the  center  in  the  case  of  a  westward  moving  tropical  cyclone.  In  this  study,  we  present 
case  studies  of  eastward  propagating  middle  latitude  (extratropical)  storms  -  in  which  the  wake 
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region  is  normally  to  the  west  of  the  center.  The  wake  region  microbarom  generation  theory  is 
discussed  in  Pnomaryov  et  al.  (1998)  and  converging  wave  trains  of  opposite  directions  were 
also  observed  in  the  wake  of  Hurricane  Bonnie  in  the  Northwest  Atlantic  Ocean  (Wright  et  al., 
1998). 

In  the  case  of  the  eastward  moving  middle  latitude  storm,  the  south  side  of  the  storm  will  contain 
the  strongest  fetch  lengths,  durations,  and  intensities  and  thus  higher  wave  heights  and  periods 
can  be  expected  on  to  the  south  of  the  center.  It  is  rare  for  an  eastward  moving  storm  to  produce 
long  wave  periods  of  1 3  seconds  or  more  on  the  north  side  of  the  center  unless  the  storm  is  a 
slow  mover,  large  in  diameter,  and/or  contains  very  high  wind  speeds.  Furthermore,  since  long 
period  swells  exhibit  higher  group  velocities  than  short  period  swells  (group  velocity  of  swells 
with  20  sec  ~  15  m/s  ~  29  knots;  5  sec  ~  3.75  m/s  ~  7  knots),  the  long  period  energy  created 
ahead  of  a  low  center  will  often  be  dispersed  before  it  has  time  to  interact  with  any  longer  period 
energy  that  may  be  created  behind  the  low  center.  This  is  especially  true  when  the  storm  is 
moving  slower  than  the  group  velocities  of  the  long  period  ocean  waves  it  produces.  Steep,  short 
period  swells  (<  8  seconds)  often  dissipate  within  a  few  wavelengths  due  to  whitecapping  and 
angular  spreading  effects.  This  can  prevent  opposing  wave  trains  of  short  periods  from 
interacting.  On  the  other  hand,  opposing  wave  trains  of  medium  periods  (8-12  seconds)  are  very 
likely  in  wake  of  both  middle  latitude  and  tropical  cyclones,  and  thus  the  infrasonic  peak  at  0.2 
Hz  is  commonly  observed  (Fig.  3).  If  there  are  long  period  swells  created  on  the  north  side  of 
the  center  (in  the  case  of  an  eastward  moving  storm),  the  amplitudes  are  normally  small  and  thus 
the  resulting  infrasonic  wave  that  results  will  also  have  small  amplitude.  Similarly,  dissipating 
shorter  period  swells  that  may  interact  will  also  produce  smaller  amplitude  infrasonic  waves. 
This  also  helps  explain  the  common  infrasonic  spectral  peak  at  0.2  Hz. 

In  this  project  the  microbarom  generation  model  of  Arendt  and  Fritts  (2000)  and  global  wave 
spectra  from  the  third-generation  Wavewatch  III  (WW3)  ocean  wave  model  of  Tolman  (1999) 
are  used  to  predict  acoustic  source  pressure  spectra.  This  will  provide  a  global  estimate  of 
probable  microbarom  generation  regions  in  the  Pacific  for  the  case  studies  presented  (early 
January  and  late  February,  2003).  The  source  modeling  results  will  also  show  the  relationship 
between  microbarom  generation  regions  and  marine  storms  and  dominant  ocean  wave 
parameters  (significant  wave  height,  peak  period,  and  mean  propagation  directions). 
Microbaroms  observed  at  the  I59US  and  other  near-Pacific  infrasound  arrays  will  be  used  to 
confirm  our  modeling  results.  The  observed  and  modeled  microbaroms  will  be  used  together  to 
help  present  a  conceptual  model  of  microbarom  generation  in  the  Pacific.  This  will  help  in  the 
characterization  of  microbarom  source  regions  observed  at  I59US  during  2003.  Finally, 
microbarom  signals  observed  in  Hawaii  during  2003  are  also  compared  to  prevailing 
atmospheric  winds  to  help  further  explain  some  of  the  seasonal  and  daily  variability  of  the 
arrivals. 
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SECTION  2 


DATA,  METHODS,  AND  INSTRUMENTATION 

The  I59US  infrasound  array  consists  of  four  Chaparral-5  differential  pressure 
microphones  with  a  frequency  range  of  0.04  to  8  Hz.  Three  of  the  sensors  are  arranged  as  a 
triangle  with  a  2  km  baseline,  with  the  fourth  sensor  near  the  center  of  the  triangle  (Fig.  2). 
Sensor  data  are  recorded  at  20  samples  per  second  by  24-bit  digitizers  and  sent  in  real  time  via 
radio  telemetry  to  the  Inffasound  Laboratory  in  Keahole  Point  on  the  west  coast  of  Hawaii.  The 
I59US  station  has  very  low  ambient  noise  levels  and  is  one  of  the  most  sensitive  stations  of  the 
IMS  because  of  its  location  in  a  dense  tropical  forest  leeward  of  Hawaii’s  massive  volcanoes. 
Porous  hose  filters  are  also  used  for  wind  noise  reduction. 

The  Progressive  Multi-Channel  Correlation  (PMCC)  algorithm  of  Cansi  et  al.  (1995)  is  the 
primary  signal  detection  system  used  at  I59US.  PMCC  uses  the  correlation  between  various 
groupings  of  three  sensors,  i,  j,  k  to  estimate  the  consistency  of  the  lag-closure  relation 

rm  =  Aty  +AtJk  +Atkl,  (1) 

where  Aty  is  the  time  delay  between  the  arrival  of  a  signal  at  sensors  /  and  j  (Cansi  and  Klinger, 
1997).  For  an  ideal  wavefront,  this  consistency  is  equal  to  zero.  A  microbarom  detection  is  only 
registered  if  the  consistency  is  below  0.5  sec  within  the  0.1 -0.5  Hz  passband  (Garces  and  Hetzer, 
2002,  2003).  The  microbarom  feature  extraction  process  uses  overlapping  windows  of  data  with 
a  length  of  90  sec  and  an  overlap  of  20  sec.  In  essence,  PMCC  is  used  to  detect  coherent 
infrasonic  energy  across  the  array,  which  allows  the  speed,  arrival  azimuth,  and  amplitude  of  the 
detected  arrivals  to  be  extracted.  Microbarom  arrival  azimuth  and  amplitude  is  the  main 
parameter  we  use  to  compare  with  storm  and  ocean  wave  characteristics  as  well  as  source 
modeling  results.  The  arrival  azimuth  of  an  inffasound  signal  is  obtained  by  determining  time 
shifts  that  yield  the  highest  cross-correlation  between  waveforms  of  the  array  elements.  Signals 
with  low  consistency,  as  well  as  a  steady  azimuth  and  trace  velocity,  are  referred  to  as  coherent 
arrivals.  Due  to  a  160  m  sensor  elevation  difference  (Fig.  2),  the  arrival  azimuth  and  apparent 
horizontal  phase  velocity  has  to  be  corrected  for  energy  incident  at  steep  angles  from  the 
horizontal  (D.  Brown,  personal  communication,  2001).  However,  most  microbaroms  arrive  near 
the  horizontal,  so  no  substantial  correction  is  needed. 

Infrasonic  power  spectral  densities  are  used  to  distinguish  peaks  in  the  microbarom  frequency 
range.  The  power  spectral  densities  were  computed  using  the  modified  periodogram  method 
(Madisetti  and  Williams,  1998)  with  a  102.4  s  Hanning  window  (211  samples)  and  a  50% 
overlap.  Power  spectra  include  the  combination  of  both  coherent  and  non-coherent  microbaroms 
at  each  frequency. 

The  WW3  ocean  wave  model  (vl.18)  is  driven  by  NOGAPS  10  m  surface  winds  and  includes 
global  ice  concentration  values  and  is  used  to  produce  realistic  ocean  wave  spectra  values  on  a 
global  1°  grid.  The  wind  and  ice  input  files  are  provided  by  the  Master  Environmental  Library 
Homepage  (MEL)  at  http://mel.dmso.mil/.  For  both  case  studies  discussed  in  this  project,  WW3 


4 


was  initialized  at  least  6  days  ahead  of  chosen  events  to  produce  an  accurate  background  ocean 
wave  field. 

The  WW3  model  outputs  wave  energy  densities,  F,  as  a  function  of  frequency,  f  and 
propagation  direction,  6.  The  wave  energy  densities  can  be  integrated  over  angle  and  frequency 
to  provide  the  total  wave  energy  E, 

E  =  [n  de[dfF{f,e).  (2) 

The  significant  wave  height  is  defined  as 

Hs  =  4 -JE  .  (3) 

WW3  outputs  wave  energy  densities  in  24  directional  and  25  frequency  bins  to  produce  600 
values  at  each  grid  point.  These  wave  spectra  are  then  used  to  calculate  acoustic  source  pressure 
by  summing  the  products  of  directly  opposing  wave  trains  at  each  frequency.  The  following 
algorithm  based  on  the  Arendt  and  Fritts  (2000)  model  was  used  to  predict  acoustic  source 
pressure  fields  on  a  global  1°  grid, 
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where  Pk  (Pa-m  )  is  the  acoustic  source  pressure  wavenumber  spectrum,  com  is  the  microbarom 
frequency  which  is  assumed  to  be  twice  that  of  the  ocean  surface  waves,  p,  c,  and  g  are  assumed 
to  be  1  kg-m'3,  340  m-s'2  and  9.81  m-s'2,  respectively.  Qm  contains  the  product  of  WW3  produced 
wave  energy  densities  traveling  in  opposite  directions  with  frequency  held  constant,  integrated 
over  all  24  directional  bins.  A  detailed  derivation  of  this  algorithm  is  given  in  Garces  et  al. 
(2003)  and  Appendix  A. 

WW3  is  considered  by  academia,  military,  and  ocean  enthusiasts  alike  to  be  one  of  the  best 
ocean  wave  models  of  its  kind  to  date.  Global  WW3  output  validates  well  when  compared  to 
both  buoy  and  satellite  altimeter  observations  (Tolman,  2002).  Tolman  (2002)  also  showed  the 
superiority  the  WW3  model  has  to  its  predecessor,  the  Wave  Action  Model  (WAM).  The  WAM 
model  was  shown  to  perform  poorly  in  large  ocean  basins  (i.e.  the  Pacific)  due  to  wave 
generation  and  propagation  errors,  while  WW3  has  made  improvements  in  these  schemes 
(Rogers  and  O'Reilly,  2002). 

Observational  and  source  modeling  results  are  first  shown  for  January  4,  2003  at  which  time  an 
intense  (sub  952mb)  mid  latitude  cyclone  was  moving  eastward  just  north  of  the  Hawaiian 
Islands,  along  with  6  other  surface  lows  scattered  throughout  the  Pacific  Basin.  The  second  case 
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study  is  from  February  22,  2003  when  the  North  Pacific  was  also  very  active  with  storm  and  high 
wave  activity.  Both  of  these  case  studies  were  chosen  primarily  because  of  the  strength  and 
coherence  of  the  microbarom  observations  during  the  events.  Microbarom  arrival  azimuths  from 
7  other  infrasound  stations  were  obtained  from  the  February  22  event,  which  allowed  for  the 
triangulation  of  the  microbarom  source  region.  Five  of  these  stations  are  part  of  the  IMS:  IS34 
(Mongolia,  China),  IS53  (Fairbanks,  Alaska),  IS57  (Pinion  Flat,  California),  IS  10  (Manitoba, 
Canada),  and  the  Hawaii  IS59  array.  Non-IMS  sites  PDIAR  (Pinedale,  Wyoming)  and  NVIAR 
(Mina,  Nevada)  also  recorded  microbaroms  from  the  North  Pacific  on  February  22. 

Surface  pressure  charts  were  generated  using  data  supplied  by  the  NCEP/NCAR  Reanalysis 
project  available  in  6  hourly  intervals  at  a  resolution  of  2.5°  (web  site  at  www.cdc.noaa.govA. 
The  NCEP/NCAR  data  set  provides  a  good  opportunity  to  examine  the  synoptic  scale  evolution 
of  storm  systems.  Surface  weather  and  wave  charts  are  compared  to  microbarom  observations 
and  source  modeling  results  to  distinguish  a  relationship  between  infrasonic  source  locations  and 
marine  storm  tracks. 

Since  meteorology  determines  how  sound  travels,  this  research  also  relates  seasonal  and  daily 
microbarom  arrivals  to  predominant  wind  directions  from  the  troposphere  to  the  mesosphere 
(Chapter  4).  In  order  to  produce  atmospheric  specifications  for  infrasonic  propagation  studies, 
the  Naval  Research  Laboratory  (NRL)  Ground  to  Space  (G2S)  model  (Drob  et  ah,  2003)  was  run 
to  produce  a  self-consistent  dataset  extending  from  January  1 ,  2003  to  March  29,  2004.  Global 
spectral  coefficients  of  wind,  temperature,  and  density  were  produced  at  6-hour  time  intervals. 
These  coefficients  have  a  triangularly  truncated  spectral  order  of  72  (T-72)  resulting  in  an 
effective  output  resolution  of  approximately  2.25  degrees.  The  G2S  winds  from  the  Earth’s 
surface  to  10  mb  (0  to  35  km)  are  based  on  analysis  output  from  the  lxl  degree  global  aviation 
model  (AVN)  from  the  NOAA  National  Center  for  Environmental  Prediction  (NCEP).  From  1 
to  .4  mb  (20  to  55  km)  the  G2S  winds  are  based  on  the  1.0x1. 5  global  assimilative  analysis  from 
the  NASA  Goddard  Space  Flight  Center,  Global  Modeling  and  Assimilation  Office  (GSFC- 
GMAO).  The  upper  atmospheric  conditions  (50-170  km)  are  specified  by  the  HWM/MSIS 
models. 
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SECTION  3 


MICROBAROM  OBSERVATIONS  AND  SOURCE  MODELING  RESULTS  FROM 
TWO  CASE  STUDIES  DURING  THE  WINTER  OF  2003 

The  I59US  station  gets  bombarded  with  coherent  microbarom  arrivals  from  250-360° 
during  the  Boreal  Winter,  which  is  where  the  strongest  storms  and  highest  seas  tend  to  be  located 
from  around  November  through  March.  Fewer  arrivals  come  from  northeast  azimuths  during  the 
winter  despite  the  common  presence  of  strong  storms  and  high  seas  northeast  of  the  Hawaii 
array.  This  appears  to  be  a  result  of  topographic  shadowing  and  advection  of  the  predominant 
winds  (discussed  in  further  detail  in  Chapter  4).  It  is  not  uncommon  for  the  daily  variability  of 
arrival  azimuth  to  exhibit  a  clockwise  rotation  as  storm  systems  travel  from  west  to  east  in  the 
North  Pacific.  Sometimes  however,  microbarom  arrivals  do  not  appear  to  follow  storm  systems. 
The  aim  of  this  chapter  is  to  compare  microbarom  observations  to  source  modeling  results, 
surface  weather  charts,  and  dominant  wave  parameters  to  better  characterize  microbarom  source 
regions. 


3.1  MICROBAROM  OBSERVATIONS  FROM  EARLY  JANUARY  2003 

Microbarom  arrivals  from  Jan.  1-3,  2003  concentrated  between  270-330°  (Fig.  7)  correspond  to 
surface  low  pressure  areas  and  high  wave  activity  in  the  region  during  that  time.  Of  particular 
interest  are  arrivals  from  Jan.  4-6,  2003.  Peak  microbarom  signals  on  the  4th  are  from  northwest 
azimuths.  The  arrival  peak  then  rotates  clockwise  through  the  5th  when  several  arrivals  from  the 
northeast  are  evident.  Coherent  energy  from  both  west  and  northeast  directions  are  recorded  on 
the  6th.  A  middle  latitude  cyclone  with  central  pressure  <  952  mb  was  located  just  northwest  of 
Hawaii  on  the  4th  and  then  moved  northeast  of  I59US  on  the  5th  (Fig.  8),  similar  to  the 
directions  the  microbaroms  were  coming  from.  This  storm  system  produced  massive  NW  swells 
along  exposed  Hawaiian  Island  shores  on  Jan.  5-6,  2003,  with  observed  breaking  wave  heights  > 
12  m  closing  out  Waimea  Bay. 

Seven  surface  low  pressure  centers  can  be  seen  in  the  Pacific  Ocean  and  its  nearby  waters  on  Jan. 
4,  2003  at  18Z  (Fig.  9).  Each  of  these  has  an  associated  peak  in  significant  wave  height  (Fig.  10) 
that  is  located  where  pressure  gradients  are  strongest.  Four  of  these  lows  are  in  the  southern 
hemisphere  but  little  to  no  coherent  microbaroms  from  southerly  directions  was  recorded  during 
this  time  frame.  These  observations  suggest  that  the  microbaroms  from  early  January  2003  are 
coming  from  North  Pacific  storm  systems  and/or  their  associated  regions  of  converging  high 
seas.  Comparison  of  microbarom  arrivals  at  I59US  during  this  time  to  dominant  wave  period 
and  direction  charts  (Fig.  11)  yielded  no  clear  relationship.  However,  it  is  nearly  impossible  to 
pinpoint  an  exact  source  location  region  from  the  arrival  azimuth  data  from  the  I59US  array 
alone.  Data  from  other  arrays  was  investigated,  but  no  similar  source  region  was  noted  as  signals 
from  storms  adjacent  to  those  stations  dominated. 

The  infrasonic  power  spectrum  on  Jan.  4,  2003  18Z  (Fig.  3)  shows  a  peak  of  0.1  Pa2/Hz  near  0.2 
Hz.  A  smaller  peak  of  .01  Pa2/Hz  is  seen  near  .135  Hz.  A  minimum  of  less  than  .001  Pa2/Hz  is 
evident  near  0.1  Hz.  The  microbarom  peak  at  0.2  Hz  is  typical  year  round  at  I59US.  A 
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microbarom  containing  a  frequency  of  0.2  Hz  would  theoretically  originate  from  two  swell  trains 
containing  periods  of  10  seconds  (0.1  Hz)  traveling  in  nearly  equal  but  opposite  directions.  Of 
particular  interest  is  the  spectral  increase  in  energy  from  0.1  to  0.2  Hz  from  Jan.  4  (Fig.  3)  to  Jan. 
6  (Fig.  12).  This  increase  corresponds  to  an  increase  in  interactions  of  ocean  waves  containing 
longer  periods  (10  to  20  sec). 

3.2  SOURCE  MODELING  RESULTS  FROM  JANUARY  4-6,  2003 

Six  out  of  seven  surface  low  pressure  centers  in  or  near  the  Pacific  Ocean  on  01/04/1 8Z  exhibit  a 
peak  in  acoustic  source  pressure  at  0.197  Hz  just  west  or  southwest  of  the  center  of  circulation 
(Fig.  13).  Lows  with  a  peak  in  acoustic  pressure  to  the  southwest  of  the  center  were  propagating 
to  the  northeast  while  lows  with  a  peak  to  the  west  of  the  center  were  propagating  more  towards 
the  east.  The  exception  is  the  surface  low  centered  in  the  extreme  NW  Pacific  (just  north  of 
Hokkaido,  Japan)  where  a  maximum  in  acoustic  pressure  is  noted  to  the  south  of  the  low  center. 
This  storm  was  propagating  just  east  of  due  north.  The  acoustic  source  pressure  peak  seemingly 
associated  with  this  low  also  appears  to  be  much  larger  and  broader  than  the  other  peaks.  Several 
other  high  source  pressure  regions  0.197  Hz  are  noted  throughout  the  Pacific  at  this  time  in 
regions  near  anticyclones  and  weak  surface  pressure  gradients. 

Another  interesting  finding  is  the  minima  in  source  pressure  at  0.197  Hz  found  in  the  immediate 
vicinity  and  just  southwest  of  the  strong  cyclone  just  NNW  of  the  Hawaiian  Islands.  The 
acoustic  peak  at  0. 1 97  Hz  is  further  upstream  of  the  center  in  this  storm  than  the  others.  Acoustic 
source  pressure  calculations  at  lower  frequencies  show  a  peak  closer  to  this  low  center  (Fig.  14). 
The  microbarom  source  at  0.122  Hz  shows  a  maximum  just  southwest  of  the  center  of  the  storm. 
Microbaroms  containing  frequencies  of  0.122  Hz  correspond  to  opposing  ocean  waves  colliding 
with  frequency  .061  Hz  (16.7  second  periods).  These  longer  wave  periods  are  a  direct 
relationship  to  the  strong  storm’s  high  winds  and  long  fetch  lengths.  Although  the  observed 
microbarom  spectral  peak  is  normally  ~  0.2  Hz,  microbaroms  of  lower  frequencies  can  be 
generated  in  cases  such  as  this  and  thus  must  be  considered. 

Modeled  infrasonic  source  pressures  with  frequency  0.135  Hz  from  Jan.  4  00Z  to  Jan  6  00Z  were 
also  inspected  (Fig.  15),  due  to  the  observed  increase  in  infrasonic  power  at  this  frequency 
during  the  time.  Of  the  3  lows  present  in  the  North  Pacific,  each  exhibited  a  wake  region  peak  in 
modeled  acoustic  pressure  at  0.135  Hz.  However,  there  were  also  other  high  acoustic  source 
pressure  regions  noted  in  the  North  Pacific  during  this  time  (some  >1000  km  from  storm 
location). 

As  the  strong  storm  just  northwest  of  the  Hawaiian  Islands  (Labeled  “L2”  in  Fig.  15)  intensified, 
increased  in  diameter,  and  pushed  towards  the  ENE/NE,  the  peak  in  acoustic  pressure  in  its  wake 
(at  0.135  Hz)  also  increased  in  amplitude  and  diameter  and  moved  ENE/NE.  The  microbarom 
signals  observed  at  I59US  during  this  time  show  a  similar  trend.  However,  not  all  of  the  arrivals 
during  this  time  match  up  with  the  wake  region  of  the  strong  low  L2.  There  are  also  arrivals 
coming  from  W  and  NW  directions  where  WW3  also  predicted  high  acoustic  source  pressure 
regions  after  low  L2  moved  N  and  NE  of  I59US. 
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The  clockwise  rotation  of  the  observed  arrival  azimuths  from  Jan.  4-6  track  well  with  the 
modeled  microbarom  sources  at  0.135  Hz.  The  other  arrivals  during  this  early  January  case  not 
associated  with  storm  wake  regions  validate  well  to  the  low  frequency  source  modeling  results  as 
well.  There  are  also  similarities  in  the  modeled  sources  at  0.2  Hz  with  the  observations,  but 
overall  the  observations  tend  to  track  better  with  the  lower  frequency  sources  in  this  case.  We 
postulate  that,  for  this  case  study,  the  low-frequency  microbarom  source  dominated  amongst  all 
possible  arrivals  because  of  its  high  energy,  large  spatial  extent,  proximity  to  the  island,  and 
longer  period.  For  these  energetic  swells,  the  prevailing  winds  may  play  a  secondary  role  to  the 
source  strength  and  arrivals  may  reach  the  station  from  any  azimuth. 

3.3  MICROBAROM  OBSERVATIONS  FROM  LATE  FEBRUARY  2003 

In  contrast  to  the  early  January  event,  microbarom  arrival  azimuth  at  I59US  from  February  19- 
25,  2003  shows  little  directional  variation.  Most  of  the  energy  between  0.1  and  0.5  Hz  during 
this  period  came  from  310°  +/-  10°  (Fig.  16).  There  were  several  eastward  propagating  surface 
lows  in  the  North  Pacific  during  this  time  frame.  This  pattern  is  unlike  the  January  3-5  pattern, 
which  showed  a  clear  clockwise  rotation  in  the  arrival  azimuths  as  a  strong  storm  system  moved 
east.  The  continuous  signal  from  -310°  (instead  of  arrivals  following  storm  systems)  during  this 
time  suggests  that  microbarom  source  generation  regions  are  not  simply  related  to  surface  low 
(Fig.  17)  or  highest  ocean  wave  (Fig.  18)  location.  No  clear  relationship  was  noted  between 
arrival  azimuth  and  peak  wave  period  or  directions  (Fig.  19)  either. 

On  February  22,  2003  low  frequency  energy  below  0.6  Hz  that  originated  from  a  region  in  the 
North  Pacific  was  recorded  on  8  different  infrasound  stations  throughout  the  Northern 
Hemisphere  and  lasted  for  several  hours.  This  exceptional  microbarom  burst  was  much  stronger 
than  normal  microbarom  events  (Fig.  5)  and  thus  was  recorded  by  eight  infrasound  arrays:  IS  10, 
Canada,  IS34,  Mongolia,  IS53,  Alaska,  IS56,  Washington,  IS57,  California,  IS59,  Hawaii, 
PDIAR,  Wyoming,  and  NVIAR,  Nevada.  Near  the  peak  of  this  inffasound  event  at  00Z  on  the 
22nd,  three  surface  low  pressure  systems  were  evident  in  the  North  Pacific.  The  first  was  an 
intense  cyclone  just  east  of  Japan,  centered  near  40N,  168E  with  a  minimum  central  pressure 
below  960mb.  A  second,  moderately  strong  and  symmetrical  surface  low  (988mb)  was  located 
near  45N,  155W  while  a  third  but  much  weaker  closed  low  was  centered  just  north  of  the 
western  Aleutian  Islands  near  55N,  172W.  The  first  two  cyclones  were  propagating  towards  the 
ENE  while  the  3rd  was  moving  slowly  and  erratically. 

The  microbarom  back  azimuths  extended  along  great  circle  routes  from  7  of  the  8  inffasonic 
arrays  during  the  peak  of  this  event  intersect  in  a  confined  region  between  25-32N  and  168- 
170W,  which  is  not  a  region  where  a  surface  low  or  peak  in  significant  wave  height  is  noted. 
The  only  station  whose  microbarom  back  azimuth  did  not  intersect  through  this  region  was  IS  1 0. 
ISlO’s  great  circle  intersects  with  IS53  just  south  of  the  central  Aleutian  Islands,  with  IS34  just 
west-southwest  of  the  first,  strong  low  and  with  IS59  near  32N,  145E.  The  observations  of  this 
event  are  discussed  in  detail  in  Bhattacharyya  et  al.  (2003). 
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3.4  SOURCE  MODELING  RESULTS  FROM  FEBRUARY  22,  2003  00Z 


Highest  predicted  acoustic  source  pressure  values  in  the  Pacific  Basin  from  00Z  on  February  22 
are  evident  in  the  Northwest  Pacific  basin.  Much  of  the  modeled  high  acoustic  source  pressures 
appear  to  be  located  in  or  near  the  wake  regions  of  the  two  surface  lows  that  were  propagating 
east-northeast  (Fig.  20).  However,  some  of  the  higher  acoustic  pressure  values  are  evident 
between  the  two  surface  lows.  This  is  associated  with  swells  with  a  westerly  component 
produced  by  the  westernmost  storm  interacting  with  swells  with  an  easterly  component  produced 
by  the  easternmost  storm.  Nearly  all  microbarom  energy  received  at  I59US  during  this  time 
frame  came  from  the  direction  of  these  modeled  peaks  in  acoustic  source  pressure.  However,  it 
is  difficult  to  tell  exactly  which  region  of  converging  ocean  waves  the  inffasound  is  coming  from 
since  the  modeled  peaks  in  acoustic  pressure  are  of  similar  azimuth  from  the  I59US  array. 

Other  maxima  in  modeled  acoustic  source  pressure  are  noted  besides  the  major  North  Pacific 
peaks:  a  broad  maxima  east  of  Hawaii  extending  north  of  the  Equator  through  50N  between 
longitudes  120  and  150W,  a  peak  just  west  of  Peru,  and  several  other  Pacific  Basin  high  pressure 
values  scattered  throughout  the  southern  hemisphere  middle  to  higher  latitudes.  Several  of  these 
are  in  regions  of  a  dominating  anticyclone  or  weak  pressure  gradients. 

The  opposing  wave  train  theory  verifies  well  when  we  compare  the  infrasound  observed  from  the 
8  stations  during  this  event  with  the  WW3  produced  source  values.  The  great  circle  paths 
associated  with  the  arrival  azimuths  at  these  stations  do  not  intersect  over  a  peak  in  significant 
wave  height  (Fig.  18)  or  near  the  center  of  the  major  surface  lows  in  the  region  (Fig.  17). 
Instead,  the  great  circles  intersect  in  regions  of  the  WW3  produced  high  acoustic  source 
pressures  (Fig.  21). 
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SECTION  4 


COHERENT  MICROBAROM  ARRIVALS  AT  I59US  DURING  2003  AND  THEIR 
RELATIONSHIP  WITH  STORMS,  TOPOGRAPHY,  AND  PREVAILING 

ATMOSPHERIC  WINDS 

Microbarom  arrival  azimuths  at  the  I59US  site  during  2002  and  2003  show  an  annual  cycle  (Fig. 
22).  During  the  months  from  June  through  September  microbarom  arrivals  generally  come  from 
east  (55-130°)  or  south  (160-220°)  directions.  The  concentration  of  east  arrivals  is  much 
stronger  than  the  south  arrivals  during  this  time.  The  months  of  October  through  March  show  an 
abundance  of  arrivals  from  230-360°  with  a  peak  from  northwesterly  directions.  The  winter 
arrivals  also  have  higher  amplitudes  than  other  seasons.  Arrivals  during  April,  May,  late 
September  and  early  October  appear  to  arrive  from  a  variety  of  different  azimuths  with  no 
distinct  peak  noted. 

Annual  arrivals  clearly  correspond  to  dominant  storm  activity  in  the  Pacific  Basin  but  are  also 
affected  by  topographic  shadowing  and  zonal  and  meridional  winds  throughout  the  atmosphere. 
Acoustic  shadowing  by  the  Mauna  Loa  and  Hualalai  volcanoes  may  explain  the  few  arrivals 
from  ~  150°  and  25°,  respectively.  Microbarom  arrival  azimuth  also  shows  a  fairly  strong 
relationship  with  the  dominant  wind  directions  in  the  troposphere,  stratosphere,  and  mesosphere 
(Fig.  23).  Some  of  the  seasonal  trends  in  the  microbarom  observations  can  be  explained  by  the 
winds  in  the  stratosphere  and  lower  mesosphere  (50-70  km),  while  some  of  the  daily  variability 
can  be  explained  by  the  tropospheric  and  lower  stratospheric  winds  (10-20  km).  The  wintertime 
W  and  NW  arrivals  relate  well  to  the  predominant  winds  at  these  levels  in  the  atmosphere,  as  do 
the  summer  time  arrivals  from  the  east.  However,  the  summer  time  arrivals  from  the  south  do 
not  correspond  well  with  the  winds.  Figure  24  shows  a  close  up  of  Figure  23  for  the  months  of 
January  and  February  of  2003.  These  are  some  of  the  most  active  swell  months  for  the  Hawaiian 
Islands.  Some  of  the  microbarom  arrival  azimuths  in  January  and  February  track  the  tropospheric 
winds  very  well,  corresponding  to  energy  that  would  be  refracted  back  to  the  ground  at  the  10- 
20km  height  in  the  atmosphere,  where  there  is  very  little  attenuation.  Likewise,  energy  ducted 
between  the  stratopause  and  the  ground  would  suffer  very  little  attenuation  and  thereby  retain 
relatively  large  amplitude.  An  interesting  feature  is  that  the  10-20km  wind  directions  above 
Hawaii  during  the  early  January  and  late  February  cases  presented  in  Chapter  3  show  little  to  no 
relationship  with  the  microbaroms  arrival  azimuth  at  IS59.  The  marine  storms  responsible  for 
these  microbarom  events  produced  high  open  ocean  wave  heights  over  a  wide  spectrum  of 
frequencies  between  0.05  and  0.1  Hz,  and  thus  produced  high  microbarom  source  amplitudes 
between  0.1  and  0.2  Hz.  High  source  amplitudes,  coupled  with  lower  attenuation  and  larger 
correlation  length  at  lower  frequencies,  would  permit  long-range  propagation  of  this  energy  over 
large  distances  and  preferential  detection  by  an  array.  As  in  the  January  4-6  case  study, 
microbarom  signals  may  overcome  the  atmospheric  conditions  along  the  propagation  path  if  the 
source  amplitude  is  sufficiently  strong  and  contains  a  lower  frequency  (<  0.2  Hz). 

A  total  of  21459  coherent  microbarom  arrivals  reached  I59US  in  2003  (Fig.  25).  9334  of  these 
arrivals  (43%),  the  vast  majority,  came  from  270-3 1 5°.  The  second  highest  45-degree  directional 
bin,  225-270°,  contained  4211  (19.6%)  microbarom  arrivals  during  2003.  2606  arrivals  (12%) 
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came  from  45-90°,  1605  (7.5%)  from  315-360°,  1599  (7.5%)  from  180-225°,  1168  (5.4%)  from 
0-45°,  549  from  90-135°,  while  only  377  arrivals  came  from  azimuths  135-180°. 

During  the  first  quarter  of  2003  (containing  months  December  2002,  January  and  February 
2003),  98%  of  the  microbarom  arrivals  came  from  azimuths  225-360°  (Fig.  26).  The  second 
quarter  of  2003  (March,  April,  May)  shows  a  similar  directional  dependence  with  81%  of  the 
arrivals  coming  from  225-360°  (Fig.  27).  However,  stronger  signatures  from  0-90°  (13%)  and 
180-225°  (4%>)  were  evident  in  Quarter  2  than  in  Quarter  1. 

A  very  different  pattern  of  microbarom  arrival  azimuth  was  observed  during  the  boreal  summer 
months  of  Quarter  3  (Fig.  28).  70%  of  the  arrivals  during  June,  July,  and  August  came  from  0- 
135°  while  26%  came  from  135-225°.  Arrivals  at  I59US  during  the  fall  months  of  Quarter  4 
(Fig.  29)  show  a  very  similar  pattern  to  the  springtime  Quarter  2  arrivals  -  with  72%  coming 
from  225-360°,  13%  from  0-90°,  and  10%  from  180-225°. 

The  observed  arrivals  during  2003  suggest  a  relationship  between  microbarom  arrival  azimuth, 
dominant  storm  activity  in  the  Pacific  Basin,  and  tropospheric,  stratospheric,  and  mesospheric 
winds  above  Hawaii.  Boreal  wintertime  arrivals  at  I59US  come  from  west  and  northwest 
directions,  while  summer  arrivals  come  primarily  from  east  and  south  azimuths.  Arrivals  during 
the  shoulder  seasons  are  more  evenly  distributed  around  the  compass.  Microbaroms  observed  at 
I59US  come  primarily  from  regions  where  storm  and  wave  activity  are  the  strongest  throughout 
the  year,  as  suggested  in  Chapter  3,  but  are  also  strongly  affected  by  the  prevailing  winds  in  the 
troposphere,  stratosphere,  and  mesosphere.  Reflections  from  coastlines  may  produce  a  portion 
of  the  arrivals  at  I59US. 
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SECTION  5 


SUMMARY,  CONCLUSIONS,  AND  DISCUSSION 
5.1  SUMMARY  AND  CONCLUSIONS 

In  this  research,  the  microbarom  generation  theory  of  Arendt  and  Fritts  (2000)  is  used  with  ocean 
wave  spectra  provided  by  the  WW3  model  to  predict  acoustic  source  pressure  fields  generated  by 
opposing  wave  trains  in  the  open  ocean.  Standing  ocean  surface  waves  are  created  when 
opposing  wave  trains  of  similar  frequencies  and  near  opposite  directions  meet.  This  interaction 
generates  a  nonlinear  pressure  perturbation  at  the  air-sea  interface  that  travels  as  an  infrasonic 
wave  through  the  atmosphere.  In  this  study,  two  detailed  case  studies  with  modeling  results 
from  the  winter  season  of  2003  are  documented.  In  addition,  all  of  the  coherent  microbarom 
arrivals  at  the  I59US  site  during  2003  are  analyzed  and  described.  The  objectives  of  this  research 
are  to  identify  and  characterize  microbarom  source  regions  in  the  Pacific  basin  and  to  show  their 
relationship  with  storm  and  ocean  wave  activity.  The  following  conclusions  are  drawn  from  a 
careful  comparison  of  microbarom  observations  with  simulated  acoustic  source  pressures  derived 
from  the  output  of  the  Wavewatch  III  model: 

1)  The  results  from  two  detailed  case  studies  show  that  observed  microbarom  arrival 
azimuths  coincide  well  with  strong  infrasonic  source  regions  predicted  by  the  WW3 
model,  suggesting  that  the  theoretical  basis  for  the  open  ocean  generation  of  microbaroms 
is  substantially  correct.  Mean  wave  parameters  such  as  significant  wave  height,  peak 
period,  and  mean  propagation  direction  are  not  effective  in  determining  microbarom 
source  regions.  Conversely,  WW3  is  only  able  to  produce  an  accurate  depiction  of 
microbarom  source  regions  generated  by  open  ocean  wave  interactions  when  the  entire 
spectrum  is  used. 

2)  Microbaroms  are  generated  wherever  ocean  surface  wave  trains  with  opposite 
propagation  directions  and  similar  frequencies  interact.  The  strongest  microbaroms  are 
often  generated  in  the  wake  regions  of  marine  storms,  where  the  amplitude  of  the 
opposing  wave  trains  is  greatest.  In  the  two  case  studies,  every  propagating  surface  low 
exhibits  a  modeled  wake-region  peak  in  source  pressure.  However,  high  amplitude 
opposing  wave  trains  can  occur  almost  anywhere  in  the  winter  hemisphere  where 
multiple  mid-latitude  storms  may  be  evident.  Thus,  high  acoustic  source  pressure  regions 
are  often  prevalent  at  a  distance  from  the  wave  producing  winds. 

3)  Opposing  wave  trains  and  thus  microbarom  source  regions  are  very  common  and  can 
occur  nearly  anywhere  in  the  ocean  at  any  given  time.  All  source  regions  may  produce 
coherent  or  non-coherent  arrivals  at  an  array  site  but  the  strongest  and  most  coherent 
signals  will  be  determined  by  1)  amplitude  and  frequency  of  opposing  wave  trains,  2) 
proximity  of  high  acoustic  source  regions  to  array  site,  3)  mesospheric,  stratospheric  and 
tropospheric  winds,  4)  thermospheric  refraction,  and  5)  topographic  shadowing. 

4)  The  majority  of  wintertime  microbarom  arrivals  at  I59US  come  from  west  and  northwest 
directions,  while  summer  arrivals  come  primarily  from  east  and  south  azimuths.  Arrivals 
during  the  shoulder  seasons  are  more  evenly  distributed  around  the  compass.  This  annual 
cycle  is  the  result  of  the  variation  of  the  most  energetic  storm/swell  source  locations 
during  the  year  along  with  the  prevailing  wind  directions  between  10-20  km  and  50-70 
km. 
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5)  Infrasound  stations  receive  coherent  arrivals  from  the  closest  and  strongest  source  whose 
propagation  path  is  favored  by  the  atmospheric  conditions,  therefore  weaker  signals  will 
be  masked,  including  those  generated  in  wake  of  distant  or  weaker  storms. 

5.2  DISCUSSION 

Past  research  and  this  research  show  that  microbaroms  often  arrive  from  regions  of  marine  storm 
regions  and  the  associated  high  seas.  This  study  presents  a  relationship  between  marine 
storms/high  seas  with  microbarom  arrivals  but  it  is  not  a  direct  relationship.  The  proposed 
generation  theory  for  microbaroms  does  not  appear  to  be  related  to  any  atmospheric  phenomena 
associated  with  the  storm  itself,  but  rather  to  the  ocean  waves  it  creates. 

Microbaroms  observed  at  I59US  from  January  3-6,  2003  appeared  to  come  from  a  region  near 
the  center  of  a  strong  surface  low  and  the  associated  peak  in  significant  wave  height  as  the  storm 
propagated  west  to  east.  In  reality,  these  arrivals  that  appeared  to  be  coming  from  the  storm 
center  were  likely  originating  just  in  the  wake  of  the  center,  where  the  swell  spectrum  was 
confused  with  multidirectional  components  created  ahead  of  and  behind  the  surface  low  center 
(Fig.  30).  However,  it  is  nearly  impossible  to  tell  the  exact  source  locations  of  the  microbaroms 
based  on  arrival  azimuth  data  from  one  infrasound  array  only.  Data  from  a  network  of  observing 
stations  during  the  February  22  REB  event  confirmed  that  microbaroms  are  not  necessarily 
generated  in  regions  of  high  significant  wave  height  or  where  a  marine  storm  is  located.  All  of 
the  back  azimuths  from  this  event  converged  over  regions  of  high  acoustic  pressure  (computed 
by  output  from  the  WW3  model),  which  were  not  collocated  with  peaks  in  significant  wave 
height  or  surface  low  location. 

Source  modeling  results  from  both  of  the  case  studies  exhibit  peaks  in  acoustic  pressure  in  the 
wake  regions  of  every  moving  marine  storm  in  the  Pacific  Basin.  The  microbarom  relationship 
with  wake  regions  is  clear  in  the  early  January  case  since  the  microbaroms  rotate  clockwise  with 
time.  However,  the  arrivals  at  I59US  from  February  19-25  show  little  directional  variation 
dependency  associated  with  moving  storm  systems.  The  dominant  arrival  azimuth  for 
microbaroms  at  I59US  during  this  time  was  near  310°.  During  the  Northern  Hemisphere  winter, 
the  seas  in  the  northwest  Pacific  Ocean  are  often  fully  developed  due  to  the  persistence  of  strong 
mid-latitude  cyclones.  The  storm  track  during  this  portion  of  February  likely  produced  a  broad, 
nearly  continuous  region  of  high  seas  converging  in  a  region  -310°  from  I59US  that  was 
producing  stronger  microbarom  signals  than  any  of  the  other  generation  regions  in  the  Pacific. 

When  multiple  propagating  surface  lows  are  evident  in  the  ocean  simultaneously,  each  storm 
may  generate  a  strong  region  of  microbaroms  in  its  wake  or  there  may  be  strong  source  regions 
between  the  storms.  Infrasound  sites  detect  coherent  signals  from  the  closest  and  strongest  source 
whose  propagation  path  is  favored  by  the  atmospheric  conditions;  therefore,  weaker  and/or  more 
distant  sources  are  masked.  Data  from  a  network  of  infrasound  sites  are  needed  to  determine  an 
exact  microbarom  source  region.  The  network  of  operational  inftasound  stations  near  the  North 
Pacific  may  allow  microbaroms  generated  in  the  wake  of  or  between  strong  storm  systems  to  be 
traced  back  to  their  source  region. 
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The  fact  that  infrasound  sites  detect  signals  from  the  closest  and  strongest  sources  favored  by  the 
atmosphere  explains  several  of  the  other  modeled  peaks  in  acoustic  source  pressure.  The  results 
from  both  February  22  and  January  4  showed  source  peaks  in  the  wake  regions  of  southern 
hemisphere  low  pressure  systems  but  little  to  no  coherent  microbaroms  were  recorded  from  the 
south  during  either  of  the  case  studies.  Arrivals  from  southern  hemisphere  sources  are  less 
obvious  because  of  the  effect  of  atmospheric  conditions  and  attenuation  along  lengthy 
propagation  paths.  Source  modeling  also  showed  several  weaker  peaks  scattered  through  the 
Pacific  under  regions  of  high  atmospheric  pressure  at  the  surface  and  in  weak  wind  regions.  At 
any  given  time  in  the  middle  of  the  ocean  there  may  be  numerous  multidirectional  and  multi¬ 
period  swell  events,  typically  more  than  three  but  less  than  20,  passing  through  a  given  point  - 
regardless  of  storm  or  dominant  swell  activity  in  the  vicinity.  This  is  illustrated  well  by 
comparing  the  frequency  and  directional  ocean  wave  spectrum  at  a  Central  Pacific  location  (Fig. 
31)  with  surface  pressure  (Fig.  9)  and  significant  wave  heights  (Fig.  10).  The  wave  spectrum  at 
this  location  (0°  N,  153.88°  W)  shows  five  different  swell  peaks  arriving  simultaneously,  but 
weak  surface  weather  patterns  are  dominating  and  no  peak  in  significant  wave  height  is  noted. 
Opposing  wave  trains  and  microbarom  generation  are  therefore  common  at  considerable 
distances  from  regions  of  strong  wave  generation.  In  particular,  both  case  studies  exhibit  a  broad 
peak  in  modeled  source  pressure  in  a  region  east  of  the  Hawaiian  Islands  north  of  the  Equator 
through  50N  between  longitudes  120  and  150W.  The  peaks  here  seem  to  be  of  only  slightly 
weaker  magnitudes  to  those  in  the  wake  regions  of  the  strong  North  Pacific  storms.  This  is 
arguably  a  result  of  the  scale  (log)  used  to  plot  these  tiny  pressure  waves  that  rarely  exceed  a  few 
mPa.  The  microbaroms  generated  in  this  region  are  expected  to  be  weaker  than  those  in  the 
Northwest  Pacific,  and  thus  were  not  recorded  on  a  coherent  basis,  but  may  indeed  arrive  as 
incoherent  signals.  There  are  also  two  other  factors  that  could  have  prevented  their  coherent 
arrival  at  I59US:  the  strong  deep  layer  westerly  flow  and  the  shadowing  from  Mauna  Kea, 
Mauna  Loa,  and  Hualalai  volcanoes.  Weak  source  regions  may  often  produce  incoherent  energy 
at  an  array  site,  but  the  strongest  and  closest  microbarom  source  region  will  generally  determine 
the  coherent  signals. 

Coherent  microbarom  signals  at  Kona  generally  come  from  south  or  east  directions  during  the 
boreal  summer.  It  is  suggested  that  the  arrivals  from  the  south  in  the  boreal  summer  originate 
from  opposing  wave  trains  in  the  wake  regions  of  eastward  propagating  extratropical  storms  in 
the  middle  latitudes  of  the  southern  hemisphere.  It  is  suggested  that  the  east  arrivals  originate 
from  1 )  tropical  weather  in  the  East  Pacific  and  2)  trade  wind  swells  reflected  from  the  Hawaiian 
Island  chain.  Incident  trade  wind  swells  are  normally  from  50-90°  (NE  to  E)  and  are  very 
common  on  the  windward  beaches  of  the  Hawaiian  Islands  due  to  a  persistent  subtropical 
anticyclone  that  dominates  the  Pacific  waters  to  the  northeast  much  of  the  year.  Sanderson 
(1993)  discussed  that  northeast  to  east  winds  (trades)  are  present  in  the  Hawaiian  Islands  from  85 
to  95%  of  the  time  in  summer,  and  from  50  to  80%  of  the  time  in  the  winter.  Reflected  swells  by 
a  land  mass  are  a  function  of  the  incident  amplitude  and  period,  along  with  the  steepness  and 
permeability  of  the  coastline  (Shore  Protection  Manual,  1973).  The  Windward  Shores  of  Hawaii 
have  an  abundance  of  rocky,  steep  coastline  so  microbarom  generation  by  reflections  is  likely 
(Fig.  32).  Trade  wind  swells  typically  range  in  period  from  8-10  seconds.  An  incident  trade 
swell  from  the  east  with  period  of  10  seconds  interacting  with  a  reflected  trade  swell  from  the 
west  would  produce  an  infrasonic  wave  with  frequency  0.2  Hz  as  predicted  by  the  frequency 
doubling  nonlinear  interactions  theory  (Arendt  and  Fritts,  2000).  Therefore,  microbarom 


15 


generation  by  reflections  may  explain  part  of  the  common  infrasonic  spectral  peak  at  0.2  Hz 
observed  at  I59US  throughout  the  year,  especially  during  summer  months. 

Occasionally,  infrasonic  power  spectra  observed  at  I59US  exhibit  split  peaks  within  the 
microbarom  threshold  (Fig.  33).  The  peak  at  ~  0.2  Hz  is  quasi-permanent  but  an  increase  in 
lower  frequency  energy  sometimes  leads  to  a  bifurcation  in  this  peak.  This  may  be  the  result  of 
two  strong  regions  of  opposing  wave  trains  interacting.  For  example,  microbaroms  generated  by 
~10  s  trade  wind  swell  reflection  may  be  occurring  simultaneously  with  microbaroms  being 
generated  by  longer  wave  period  swells  associated  with  tropical  cyclones  in  the  East  Pacific.  The 
trade  swell  reflection  would  lead  to  a  peak  at  0.2  Hz  while  the  tropical  cyclones  produce  a  split 
peak  at  lower  frequencies.  This  situation  can  also  occur  in  the  winter  when  multiple  mid-latitude 
cyclones  in  the  North  Pacific  produce  numerous  regions  of  converging  wave  trains  of  varying 
periods  (as  shown  in  the  early  January  case).  Thus  it  is  possible  to  explain  the  structure  of  ocean- 
related  infrasound  spectra  observed  in  Hawaii  by  attributing  (1)  a  broad  2-5  Hz  energy  peak  to 
breaking  ocean  waves  (surf)  in  the  littoral  zone,  (2)  0.2  Hz  energy  to  the  predominant  10  s  swell 
energy  and  (3)  0.1 -0.2  Hz  energy  to  12-20  s  period  swells  associated  with  powerful  storms.  The 
surf  peak  depends  on  the  swell  energy  arriving  near  the  array  site,  and  in  Hawaii  surf  sound 
appears  to  be  propagated  in  the  boundary  layer  and  the  troposphere.  The  coherent  component  of 
the  main  microbarom  peak  at  0.2  Hz  depends  on  the  prevailing  winds  and  the  amount  of  swell 
energy  in  the  ocean.  However,  the  energy  of  the  main  microbarom  peak  is  the  sum  of  coherent 
and  incoherent  components,  and  the  amplitude  of  this  peak  may  be  related  to  the  diurnal  tides  in 
the  upper  atmosphere  (Rind,  1978).  If  atmospheric  winds  do  not  support  tropospheric  or 
stratospheric  waveguides,  a  low  trace  velocity  could  correspond  to  refraction  in  the  mesosphere 
and  lower  thermosphere.  Coherent  microbaroms  in  the  -0.12-0.2  frequency  range  may  arrive 
from  any  direction  where  energetic  long-period  swells  are  interacting,  and  under  certain 
circumstances  may  also  be  used  to  study  the  mesosphere  and  lower  thermosphere  (Garces  et  al., 
2002).  Returns  from  the  upper  atmosphere  would  explain  southerly  arrivals  associated  with  the 
powerful  southern  hemisphere  swells  during  the  boreal  summer. 

Microbarom  arrivals  were  shown  to  be  concentrated  from  north  and  west  directions  during  the 
Boreal  winter.  Arrivals  with  a  northerly  component  are  hypothesized  to  originate  from  non¬ 
linear  interactions  of  ocean  swells  generated  by  eastward  moving  extratropical  storms  in  the  mid¬ 
latitudes  of  the  North  Pacific.  Arrivals  from  west  or  west-southwest  directions  in  the  boreal 
winter  are  likely  from  Western  Pacific  tropical  cyclones,  including  those  that  recurve  during 
extratropical  transition  (when  tropical  cyclones  lose  their  warm  core  characteristics  while 
transitioning  into  a  cold  core  mid-latitude  cyclone).  The  west  arrivals  during  the  winter  could 
also  be  generated  by  maturing  cold  core  weather  systems  as  well.  Swell  convergence  associated 
with  the  Intertropical  Convergence  Zone  (ITCZ)  can  also  produce  microbarom  arrivals  at  I59US 
throughout  the  year  from  southeast,  south,  or  southwest  directions.  Arrivals  during  fall  and 
spring  months  were  shown  to  come  from  a  variety  of  directions  with  no  distinct  trend  noted,  due 
to  the  randomness  of  storm  systems  in  the  Pacific  Basin  during  shoulder  seasons. 

The  varying  atmospheric  conditions  throughout  the  year  are  also  expected  to  have  a  great 
influence  on  the  microbarom  arrivals.  Microbaroms  are  believed  to  travel  through  a  great  extent 
of  the  atmosphere  instead  of  directly  on  the  horizontal  (Gossard  and  Hooke,  1975).  This  leaves 
the  sound  waves  vulnerable  to  numerous  influences  from  the  varying  atmospheric  conditions 
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during  the  year.  The  change  in  temperature  with  height  will  affect  the  propagation  of  the  sound 
waves,  as  the  phase  speed  of  sound  is  proportional  to  the  square  root  of  the  temperature.  I59US 
may  be  more  inclined  to  receive  microbaroms  from  the  west  and  northwest  in  the  winter  and 
from  the  east  and  northeast  during  summer  because  of  the  advection  of  the  sound  waves  by  the 
mean  flow  as  show  in  Chapter  4.  Figures  23  and  24  suggest  that  seasonal  trends  in  microbarom 
arrival  azimuth  are  associated  with  stratospheric  wind  directions  while  daily  variability  in 
arrivals  may  be  explained  by  upper  tropospheric  wind  patterns.  In  addition,  substantial  diurnal 
variability  may  be  introduced  by  diurnal  boundary  layer  fluctuations  (Fig.  34).  However,  not  all 
of  the  arrivals  can  be  explained  solely  by  wind  direction.  The  strong  microbarom  signals 
observed  at  I59US  during  the  case  studies  presented,  for  instance,  show  little  correlation  with  the 
winds  over  Hawaii.  This  suggests  that  a  strong,  close  microbarom  source  region  may  produce 
coherent  energy  at  an  array  site  regardless  of  the  atmospheric  conditions  in  the  region.  The  issue 
of  atmospheric  wind  effects  on  annual  microbarom  arrivals  at  I59US  is  discussed  in  detail  in 
Garces  et  al.  (2004). 

5.3  CONCLUDING  REMARKS  AND  FUTURE  WORK 

Although  in  general  the  continuous  microbarom  noise  is  uncorrelated  at  I59US,  distinct 
microbarom  packets  can  be  detected  in  Hawaii.  When  designing  an  array  it  may  be  useful  to 
retain  the  coherence  of  microbaroms,  as  then  it  is  possible  to  remove  this  contribution  from  a 
signal  of  interest. 

This  research  computed  acoustic  source  pressures  calculated  by  colliding  ocean  waves  given  by 
the  global  1°  WW3  model.  There  will  certainly  be  subgrid  scale  regions  that  contain  opposing 
wave  trains  and  thus  microbarom  generation  regions.  The  next  step  of  this  project  will  involve 
transporting  the  WW3  model  to  the  Maui  High  Performance  Computing  Center  (MHPCC)  to 
refine  the  resolution  of  both  the  wave  and  microbarom  source  modeling.  A  high  resolution 
bathymetry  file  and  shallow  water,  near  shore  ocean  physics  would  need  to  be  incorporated  into 
the  ocean  modeling  to  show  the  microbarom  source  regions  associated  with  reflections  from  a 
coastline. 

Producing  an  accurate  model  of  the  microbarom  source  is  only  the  first  step  in  a  very  complex 
problem.  The  modeled  source  pressures  need  to  be  propagated  across  the  globe  with  atmospheric 
specifications  (wind,  temperature,  pressure,  and  density),  topography,  attenuation,  and 
propagation  time  taken  into  account  in  order  to  define  microbarom  noise  at  an  array  site.  Global- 
scale  infrasound  propagation  modeling  is  beyond  the  scope  of  this  report. 

The  results  presented  in  this  document  show  there  is  potential  for  using  microbaroms  to  track 
ocean  wave  activity.  Buoy  and  ship  observations  are  scarce  while  wave  parameters  from  satellite 
interpolations  are  not  continuous  because  they  are  obtained  through  satellites  on  polar  orbits. 
Incidentally,  infrasonic  energy  in  the  2-5  Hz  range  associated  with  surf  breaking  has  been  shown 
to  be  closely  correlated  to  significant  wave  heights  recorded  by  a  nearby  buoy.  This  also  is 
beyond  the  scope  of  this  paper  but  is  discussed  in  Garces  et  al  (2003).  Microbaroms  contain 
information  about  the  wave  periods  and  amplitudes  of  the  generating  open  ocean  waves.  Wave 
period  is  often  a  parameter  that  goes  overlooked  in  operational  forecasting,  but  in  fact  is  a  crucial 
component  to  the  prediction  of  breaking  wave  heights. 
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APPENDIX  A 


DERIVATION  OF  BASIC  RELATIONSHIPS  FOR  THE  SURFACE  WAVE  SPECTRUM, 
ACOUSTIC  SOURCE  FUNCTION,  AND  INFRASONIC  AMPLITUDE  SCALING  FOR 

MICROBAROM  SIGNALS 

A.l  SURFACE  WAVE  SPECTRUM 

The  Wavewatch  3  model  outputs  the  variance  density,  F,  of  the  surface  wave  field  as  a  function 
of  frequency,  /  and  propagation  direction,#.  The  variance  density  can  be  integrated  over  angle 
and  frequency  to  provide  the  total  wave  energy  E, 

E=l’de[(tfF{f,e).  (Al) 

The  significant  wave  height  is  defined  as 

H,  =  4 JE  .  (A2) 

And  thus  the  variance  density  has  units  of  m  /(rad*Hz)  and  it  is  a  measure  of  the  energy  in  the 
surface  wave  field.  The  phase  of  each  wave  component  is  assumed  to  be  random. 

Arendt  and  Fritz  (2000)  start  with  a  prescribed  amplitude  spectrum  of  the  ocean  surface  wave 
field  given  by  the  cosine  Fourier  transform 

g(x,y,t)=  £ dkx  £ dkyA(kx , ky )cos(a>t -kxx-kyy  +  </>)  (A3) 


=  £  dkx  £  dky  A(kx ,  ky  )cos [cat  -  kxx-  kyy  +  <p) 

Thus  we  wish  to  find  a  functional  relationship  between  the  variance  density  F  (also  known  as  the 
directional  spectrum)  and  the  Fourier  coefficient  A,  which  has  to  have  units  of  m3. 

From  Kinsman  (1965),  Eq.  (8.3:1) 

?(*,.>', r  cos(ftrt  -  kxx  -  kyy  +  (p\jF(o),0)d6do)  (A4) 

This  integral  representation  reminds  us  that  we  are  dealing  with  a  Gaussian  process,  and  is  the 
subject  of  an  involved  discussion  by  Kinsman.  Since  the  WW3  model  is  defined  in  a  grid  of 
frequencies  and  azimuths,  we  estimate  our  wavenumber  spectrum  as: 

g(x,y,t)«  |°  £cos (cot-kxx-kyy  +  <t>]^^^-d6dco 

The  first  transformation  only  concerns  the  change  of  variables  from  frequency  to  wavenumber, 
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y,t) «  |°  £  cos(^  -  kxx  -  kyy  +  <f>)ps  A^jp-dedk 


(A5) 


where  cg  is  the  group  velocity.  The  second  transformation  involves  a  change  of  coordinates  in 
wavenumber  space  from  polar  to  rectangular, 

g(x,y,t)~  £  dkx  £  dky  cos  (cot  -kxx-  kyy  +  J  (A6) 

For  short-period  waves,  the  phase  velocity  c*  and  group  velocity  cg  are  given  as 


„  _  v  _  \S  _  S  _  1 
'  Jfc  U  <y’  Cg  2 C,t  ’ 


(A7) 


and  the  dispersion  relation  is  co2  =  gk  .  Comparison  of  Eq.(A6)  with  Eq.(A3)  yields 

Rm  j 


(A8) 


■>'/  &  V  AOAco 

since  the  variables  in  the  radical  have  units  of  (m  s),  we  see  that  /f  has  units  of  m  ,  as  desired. 

A.2  ACOUSTIC  SPECTRUM 

From  Arendt  and  Fritz  (2000)  the  pressure  field  radiated  by  a  three  dimensional  surface  wave 
field  into  an  isotropic  half  space  may  be  expressed  as 

co  oo  oo 

jdK,  ldKyldK,P„(K„Ky,Kz)exp(aJ-K»R  +  ^,  (A9) 


-OO  —00 


where  K  is  the  acoustic  wavenumber,  com  is  the  angular  frequency  for  the  sound  field,  and  R  is 
the  spherical  radius  from  the  source  location  to  the  receiver.  The  spectrum  Pk  is  given  by 


P,  = 


pc_ 

64g: 


■KZK'  1 


(A10) 


where  k  and  k ’  correspond  to  the  second  order,  propagating  solutions  of  interacting  wave  fields 
with  near  equal  angular  frequencies  co  and  co'  and  opposite  directions,  and 


K 2  =  K]  +  K2  +  K 2 

x  y  z 


(All) 
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Kx  =  kx  +  k'x ,  E,x  =  kx  -  k'x 
Kv  =k  +k'  ,  £=k-k\ 


Kz  = 


CO  +  (o' 

C  ) 


-K-K  • 


The  highest  spectral  amplitude  will  occur  when  Kx  =  Ky  =  0,  or  k’  =  -k,  and  co=  co\  and 


P'  =  77~7 K ' - k>  M-  ■ K  W  (A  1 2) 

64  g  0 


K  =  KZ  = 


2  co 
c 


con 


c 


Z*  =  2kx 


(A  13) 


Where  #is  now  the  directional  angle  for  k  specified  in  equation  (1).  Using  Equations  (6)  and  (8), 
as  well  as  com  =  2<o, 


c  ' 


A(kx,ky)A(-kx-ky)=  -f  *„-*,)  =  -4-  ylF(kx,ky]F(-kx-ky).  (A  14) 


k  J 


V®-  j 


the  peak  spectrum  amplitude  then  becomes 


4  <y„ 


2/r 


peg 

4ar 


(A  15) 


Ql  -  =  2 JVf(/,<F(/,0 +  *))<*?  (A16) 

0  0 

the  integral  can  be  readily  evaluated  through  either  the  trapezoidal  rule  or  the  parabolic  rule.  For 
an  even  number  of  intervals  n  (odd  number  of  discrete  points  Xj,j  =  0 ,n),  the  parabolic  rule  states 

\f{x)dx  « |  \f{xQ )  +  4  f(xx )  +  2  f(x2 )+...  +  4 /(*„_, )  +  /(*„ )] , 

where  h  =  (b-a)/n.  If  n  is  an  odd  number  (even  number  of  points),  we  can  use  the  trapezoidal 
rule. 
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J/M*  *  |  [/t*o )  ■ +  2/(*. )  +  2/(x2  )+•••  +  2/(x„_, )  +  /(*. )]  ■ 

a 

After  this  integration  is  complete,  we  have  completed  the  transformation  from  ocean  to  sound 
waves  and  it  is  practical  to  use  the  acoustic  microbarom  frequency  com  =  2 co . 

A.3  PRESSURE  SCALING  WITH  RANGE,  INCLUDING  THE  EFFECTS  OF  THE 
STATOSPHERIC  WINDS 

The  atmosphere  is  a  complex  environment,  and  for  the  purpose  of  microbarom  amplitude 
determinations  it  behooves  us  to  start  with  a  simple  amplitude  scaling  relationship  that  will  allow 
us  to  estimate  relative  acoustic  levels  at  a  receiver.  Stevens  et  al.  (2002)  provide  an  evaluation  of 
various  scaling  laws  that  have  been  derived  from  nuclear  tests,  and  we  use  the  expression  that 
provided  the  best  fit: 

P  =  -^r  100019v 

R]i6 

where  the  reference  pressure  Prej  is  relative  to  the  source  location,  the  range  R  is  in  km,  and  the 
stratospheric  wind  speed  v  at  50  km  along  the  propagation  direction  is  given  in  m/s.  A  source 
magnitude  may  be  defined  as 


Mi  -  logio(P ref), 


and  used  to  plot  the  microbarom  source  spectrum  amplitude  as  a  function  of  global  location. 
However,  this  expression  does  not  account  for  the  substantial  travel  time  it  may  take  for  the 
signal  to  reach  the  station,  and  the  fact  that  the  ocean  wave  field  would  have  changed  during  that 
travel  time.  Although  the  Stevens  formula  was  used  for  preliminary  evaluations  of  the  predicted 
infrasound  field,  the  results  were  not  satisfactory  and  are  not  presented  in  this  document. 
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♦  IMS  Stations 

♦  US  IMS  Stalioos 


Figure  A-l.  Infrasound  network  of  the  International  Monitoring  System. 
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Figure  A-2.  Geographic  location  of  infrasound  array  near  Kona,  Hawaii  (19.5915,  155.8936 
W).  Blue  dots  on  bottom  figure  represent  microphone  location. 
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Power  Spectral  Density,  Jan  04  2003,  08:00 


Frequency  (Hz) 

Figure  A-3.  Acoustic  power  spectral  density  observed  at  I59US  site  on  January  4,  2003  18Z 
(8am  Local).  Colored  lines  represent  power  observed  by  the  four  components  of  the 
infrasonic  array.  Power  (PaVHz)  is  on  the  vertical  axis,  acoustic  frequency  (Hz)  is  on  the 
horizontal  axis.  Blue  shaded  region  represents  microbarom  frequency  range.  The  ~0.2  Hz 
peak  is  common  throughout  the  year  at  I59US,  which  theoretically  corresponds  to  ocean 
waves  of  10  sec. 
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Figure  A-4.  Graphical  PMCC  window  from  Jan.  4,  2003  microbarom  event. 
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Figure  A-5.  Graphical  PMCC  window  from  Feb.  22, 2003  microbarom  event. 


Figure  A-6.  Simplified,  idealized  flow  around  tropical  cyclone  (or  mid-latitude  low)  “Z” 
moving  west  at  20kt  in  northern  hemisphere.  In  the  top  figure  (Time  1),  notice  that  Fetch  X 
is  pointed  in  the  direction  of  point  A.  In  the  bottom  figure  Fetch  Y  is  pointed  in  the 
direction  of  point  A  as  well,  but  a  different  Time  2.  Note  that  Fetch  X  and  Y  are  nearly 
identical  in  magnitude  and  fetch  length  and  also  pointed  in  nearly  opposite  directions. 
Therefore,  the  ocean  waves  these  two  fetches  create  will  contain  nearly  identical  amplitudes 
and  frequencies  and  will  propagate  in  nearly  equal  but  opposite  directions.  Amplitude  will 
decay  but  frequency  and  direction  will  remain  constant. 
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Time  3 


Tropical  Cyclone  ”2"  still  propagating  W 
20kt.  Swells  created  at  earlier  times  from  2 
are  beginning  to  interact  at  point  A.  Note 
that  Fetch  X  (time  1)  and  Fetch  Y  (time  2) 
were  nearly  equal  but  In  opposite  directions. 


SwcH*  from  Fetch  X 
creeled  at  Time  1  arriving 
ot  point  A  at  Tim©  3 


Swells  from  Pole Fi  Y 
crested  si  Tims  2  arriving 
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Figure  A-7.  Swell  trains  of  nearly  identical  frequencies  and  in  nearly  equal  but  opposite 
direction  meet  at  point  A  during  Time  3.  This  scenario  supports  efficient  generation  of 
Infrasound  in  the  0.1-0.5  Hz  band.  Regardless  of  the  amplitudes  of  the  interfering  wave 
trains,  acoustic  radiation  will  still  occur  (Arendt  and  Fritts,  2000). 
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Figure  A-8.  Arrival  azimuth  of  coherent  microbaroms  at  I59US  from  January  1-7,  2003. 
Arrival  azimuth  is  on  the  vertical  axis  going  clockwise  from  north  (0°  to  360°). 
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Figure  A-9.  Analyses  of  surface  pressure  (mb)  and  surface  winds  (kt)  from  January  3, 
2003  18Z  through  January  5,  2003  06Z. 
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Figure  A-10.  Analysis  of  surface  pressure  (mb)  and  isotachs  (kt,  shaded)  on  January  4, 
2003  18Z. 
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Figure  A-ll.  WW3  Significant  wave  heights  (m,  shaded)  and  mean  propagation  direction 
vectors  (towards)  for  January  4,  2003  18Z.  Black  star  represents  location  of  I59US  array. 
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Figure  A-12.  WW3  Peak  periods  in  seconds  (shaded)  and  directions  for  January  4,  2003 
18Z. 
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Power  Spectral  Density,  Jan  06  2003,  08:00 


Figure  A-13.  Acoustic  power  spectral  density  observed  at  I59US  site  on  January  6,  2003 
18Z  (8am  Local).  Colored  lines  represent  power  observed  by  the  four  components  of  the 
infrasonic  array. 
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Figure  A-14.  Base  10  logarithm  of  the  magnitude  of  the  acoustic  source  pressure  spectrum 
(Pa  *  m3)  with  frequency  0.197  Hz,  corresponding  to  ocean  waves  interacting  with  periods 
of  ~  10  s  (produced  by  equation  4  from  spectral  output  given  by  WW3)  on  01/04/03  18Z. 
Location  of  surface  low-pressure  centers  indicated  with  “L”  and  storm  propagation 
directions  are  indicated  with  arrows. 
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Figure  A-15.  Base  10  logarithm  of  the  magnitude  of  acoustic  source  pressure  spectrum  (Pa 
*  ms)  with  frequency  0.122  Hz,  corresponding  to  ocean  waves  interacting  with  periods  of  ~ 
16.4  s  on  01/04/03  18Z.  Center  of  strong  storm  north  of  Hawaii  indicated  by  “L”. 
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Figure  A-16.  Acoustic  source  pressure  (Pa  *  m3,  log  base  10)  with  frequency  0.135  Hz 
from  Jan.  4-6,  2003.  Corresponds  to  ocean  waves  interacting  with  periods  of  ~  14.8  s.  The 
three  lows  in  the  North  Pacific  during  this  time  are  tagged  Lj,  L2,  L3. 
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Figure  A- 17.  Arrival  azimuth  of  coherent  mierobaroms  at  IS59  from  February  19-25, 
2003.  Arrival  azimuth  is  on  the  vertical  axis  going  clockwise  from  north  (0°  to  360°). 
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Figure  A-18.  Analysis  of  surface  pressure  (mb)  and  isotachs  (kt,  shaded)  on  February  22, 
2003  00Z. 
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Figure  A-19.  WW3  Significant  wave  heights  (m,  shaded)  and  mean  propagation  direction 
vectors  (towards)  for  February,  22  2003  00Z.  Black  star  represents  location  of  I59US  array. 


A-23 


Figure  A-20.  WW3  Peak  periods  (shaded)  in  seconds  and  directions  for  February  22,  2003 
00Z. 
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Figure  A-21.  Base  10  logarithm  of  the  magnitude  of  the  acoustic  source  pressure  (Pa  *  m3) 
with  frequency  0.197  Hz,  corresponding  to  ocean  waves  interacting  with  periods  of  ~  10  s 
(produced  by  equation  4  from  spectral  output  given  by  WW3)  on  02/22/03  00Z.  North 
Pacific  low  pressure  centers  and  propagation  direction  indicated  with  L’s  and  arrows, 
respectively.  The  low  just  north  of  the  Aleutian  Islands  was  moving  very  little. 
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Figure  A-22.  As  in  Fig.  20  but  reshaded  to  clarify  peak  regions.  Red  shading  represents 
high  pressure  values  while  blue  shading  represents  low  pressure  values.  Great  circle  paths 
overlaid  correspond  to  the  measured  microbarom  azimuths  of  Feb.  22,  2003  00Z  REB 
event. 
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Figure  A-23.  Time  series  of  coherent  microbarom  arrivals  at  the  I59US  site  for  2002  and 
2003.  Volcano  shadow  zones  overlaid  with  dotted  lines.  Note  the  decrease  in  amplitude 
observed  during  Summer. 
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Figure  A-24.  Coherent  microbarom  arrival  (black  circles)  and  wind  arrival  azimuths, 
clockwise  from  North,  at  the  Hawaii  array.  The  transparent  circles  with  the  red  rim 
represent  the  winds  between  50-70  km  and  the  green  circles  represent  the  winds  between 
10-20  km  evaluated  at  18UT.  The  dominant  wind  directions  match  the  seasonal  variability 
for  some  of  the  arrivals,  except  for  the  arrivals  from  the  Southern  hemisphere  during  the 
Austral  winter.  These  S  swells  are  large,  consistent,  and  powerful,  and  may  overwhelm  the 
10s  period  swell  energy. 
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Figure  A-25.  Same  as  Figure  A-5  but  only  for  January  and  February,  and  showing  a 
match  between  the  tropospheric  wind  direction  and  a  number  of  the  microbarom  arrival 
azimuths. 
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Figure  A-26.  Total  microbarom  arrivals  at  I59US  during  2003.  Arrival  azimuth  in  degrees 
is  located  on  the  x-axis,  total  arrivals  on  the  y-axis. 
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Figure  A-27.  Total  microbarom  arrivals  at  I59US  during  peak  boreal  winter  months 
(December  2002,  January  &  February  2003).  Arrival  azimuth  in  degrees  is  located  on  the 
x-axis,  total  arrivals  on  the  y-axis. 
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Figure  A-28.  As  in  Figure  A-5  for  months  March,  April,  May  2003. 
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Figure  A-29.  As  in  Figure  A-5  for  months  June,  July,  August  2003. 
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Figure  A-30.  As  in  Figure  A-5  for  months  September,  October,  November  2003. 
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Figure  A-31.  Frequency,  directional  ocean  wave  spectrum  for  a  grid  point  (38.00N, 
170.00W)  at  18Z  on  Jan.  4,  2003  in  the  wake  region  of  the  strong  marine  storm  shown  in 
Fig.  5.  Frequency  (Hz)  decreases  towards  the  center,  wave  energy  scale  (m2/Hz*Deg)  on  the 
right  hand  side. 
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Figure  A-32.  As  in  Figure  24  but  for  Central  Pacific  location  0.00N,  153.88W  on  January 
4,  2003  18Z  (region  of  weak  surface  winds). 
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Figure  A-33.  Reflected  and  incident  trade  wind  swells  creating  standing  wave  along  a  steep 
coastline  of  the  Hawaiian  Islands. 
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Figure  A-34.  Acoustic  power  spectral  density  observed  at  I59US  site  on  June  16,  2003  18Z 
(8am  Local).  Colored  lines  represent  power  observed  by  the  four  components  of  the 
infrasonic  array. 
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Figure  A-35.  Spectrogram  for  all  of  2003  (upper  panel)  and  January  2003  (lower  panel)  at 
IS59  in  the  0.1-1  Hz  frequency  band.  The  upper  panel  shows  the  annual  variability  in  the 
microbarom  amplitude  and  frequency,  whereas  the  lower  panel  shows  the  diurnal 
variability  introduced  by  boundary  layer  wind  fluctuations. 
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